Monogamy of quantum correlations reveals frustration in quantum Ising spin system: 

Experimental demonstration 
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We report a nuclear magnetic resonance experiment, which simulates the quantum transverse 
Ising spin system in a triangular lattice and further show that the monogamy of quantum corre- 
lations can be used to distinguish between the frustrated and non-frustrated phases in the ground 
state of this system. Adiabatic state preparation methods are used to prepare the ground states of 
the spin system. We employ two different multipartite quantum correlation measures to analyze the 
experimental ground state of the system in both the frustrated and non-frustrated regimes. In par- 
ticular, we use multipartite quantum correlation measures generated by monogamy considerations of 
negativity, a bipartite entanglement measure, and that of quantum discord, an information-theoretic 
quantum correlation measure. As expected from theoretical predictions, the experimental data con- 
firms that the non-frustrated regime shows higher quantum correlations compared to the frustrated 
one. 



I. INTRODUCTION 

The last two decades have witnessed a phenomenal 
growth in the domain of employing quantum correlations 
[H |2] in quantum information processing tasks, with ap- 
plications ranging from quantum communication [5] to 
measurement-based quantum computation |4 and deter- 
ministic quantum computation with a single qubit [5]. 
Such advantages of using quantum devices over their 
classical counterparts have been experimentally demon- 
strated in several physical systems, which include ion 
traps [6], photonic devices [7], nuclear magnetic reso- 
nance (NMR) [8] , and cold atoms in optical lattices [9] . In 
particular, the Shor's factorization algorithm [10] . which 
provides a polynomial time quantum algorithm for find- 
ing prime factors of integers, with all known classical 
algorithms for the same task being non-polynomial time 
ones, has revolutionized the developments in the fields of 
quantum computation and quantum cryptography. The 
first experimental demonstration of the Shor's algorithm 
was performed by NMR techniques 

Recently, quantum correlations have also been em- 
ployed beyond its traditional ambit of quantum compu- 
tation and quantum information. In particular, quan- 
tum correlation has been proposed to be a universal de- 
tector of quantum phases in quantum many-body sys- 
tems [21 [T21 [13] . Detecting phases of many-body systems 
is an important task with fundamental as well as tech- 
nological implications [2]. Frustrated cooperative phe- 
nomena form one of the centerstages of the research in 
many-body physics. Frustration in a many-body Hamil- 
tonian appears when it is not possible to simultaneously 
minimize each of the energy bonds in the Hamiltonian 
independently, and may appear as a result of compet- 
ing interactions due to the geometry of the lattice or 
due to incommensurate values of the coupling strengths 



[15] . Frustrated systems usually possess hugely degener- 
ate ground states and a rich phase diagram ranging from 
quantum spin liquids to resonating valence-bond states 
|16j . Very recently, theoretical studies have shown that 
entanglement can be an effective tool for investigating 
frustrated quantum systems [TTJ, [TB] . Frustrated interac- 
tions had been known to be present in several solid state 
systems [16] . Recent experimental breakthroughs have 
made it possible to engineer frustrated spin models in ul- 
tracold atoms in optical lattices []jj], trapped ions |20) . 
NMR [21;, etc, and have led to the possibility of observ- 
ing the effects of entanglement in the different phases of 
frustrated spin models in the laboratory. 

In the present work, we employ multipartite quantum 
correlations to distinguish the frustrated phases from the 
non-frustrated ones. There are many ways by which mul- 
tipartite quantum correlations can be quantified [TJ [2] 
and even though their properties can vary widely, there 
are some distinct connecting themes. One of them is the 
"monogamy" of bipartite quantum correlations, which 
broadly demands that if two parties are strongly quan- 
tum correlated, they cannot have a significant amount of 
the same with a third party [22H21]- 

There is a two-fold aim of the present work. First 
of all, we want to experimentally observe the effect of 
monogamy of bipartite quantum correlations in a mul- 
tiparty quantum system. Secondly, we wish to apply 
the concept of monogamy of such correlations to distin- 
guish between frustrated and non-frustrated phases in a 
lattice of quantum spins. To attain our goal, we pre- 
pare a transverse Ising spin system |25j in a triangular 
configuration by using NMR techniques. It undergoes 
a transition from a non-frustrated phase to a frustrated 
one when its coupling strength is varied from a ferro- 
magnetic to an antiferromagnetic regime. We study the 
multipartite quantum correlations of the ground state to 
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detect the two different regimes. Specifically, we con- 
sider multiparty quantum correlation measures, gener- 
ated from monogamy studies of bipartite quantum cor- 
relations [22 -24 . In analyzing the experimentally gen- 
crated ground state, we employ monogamy of (i) neg- 
ativity [24] [26j [27], which is a bipartite entanglement 
measure, and of (ii) quantum discord [281 129| , which is 
an information-theoretic quantum correlation measure. 
For such investigations, we initially prepare the spin sys- 
tem in the ground state of the transverse field Hamil- 
tonian, which is a product state and then adiabatically 
drive it to both frustrated and non-frustrated regimes 
in such a way that the system remains in the ground 
state of the instantaneous Hamiltonian. We then calcu- 
late the multipartite correlations of the ground state by 
performing quantum state tomography at different time 
intervals. Moreover, coinciding with the theoretical sim- 
ulations, we observe that the non-frustrated regime has 
higher multipartite quantum correlations than the frus- 
trated one. The transition point from the non-frustrated 
to the frustrated regime is well indicated by the vanishing 
monogamy relation. 

The paper is arranged as follows. In Sec. [TTJ we present 
the multipartite quantum correlations that are employed 
in analyzing the ground state. We describe the frustrated 
and non-frustrated quantum Ising spin system in Sec. III. 
The experimental results are presented in Sec. |IV| and we 
conclude in Sec. [V] 

II. MULTIPARTITE QUANTUM 
CORRELATIONS 

Quantum correlations of a multipartite quantum state 
can be quantified in a variety of approaches. A promi- 
nent one among them is to use the concept of monogamy 
of bipartite quantum correlations [24l EH [29] . In a tripar- 
tite scenario, monogamy of a bipartite quantum correla- 
tion restricts the amount of that correlation which can 
be shared between the three parties. Let us suppose that 
three parties, Alice, Bob, and Charu, denoted respec- 
tively as 1, 2, and 3, share a tripartite quantum state. 
Monogamy of quantum correlations implies that if 1 has 
substantial quantum correlations with 2, it can have only 
a restricted amount of the same with 3. To state it more 
precisely, for a bipartite quantum correlation measure Q 
and a tripartite quantum state P123, one can introduce a 
quantity, known as "monogamy score for Q" and denote 
as 5q, given by 

8q{P12z) = Ql(23) - Ql2 - Ql3, (1) 

where Qi(23) is the bipartite quantum correlation mea- 
sured for the state P123 in the 1 : 23 partition. Q 12 is 
the same measure for the reduced state P12 — tr3pi23, 
and similarly for Q13. The three-party quantum state 
P123 is said to be monogamous with respect to the bipar- 
tite quantum correlation measure Q, if ^3(^123) > 0. A 
measure for which all tripartite quantum states produce 



non-negative monogamy scores is said to be "monoga- 
mous" . 

The monogamy score, and its properties, will cer- 
tainly depend on the bipartite quantum correlation 
measure that is employed to define the score. There 
are many ways in which one can conceptualize bipar- 
tite quantum correlations, and correspondingly there 
are many bipartite quantum correlation measures [TJ 
[2J. Broadly, these measures are defined within two 
paradigms, viz. the entanglement-separability paradigm 
and the information-theoretic one. We consider bi- 
partite quantum correlation measures from both these 
paradigms to define monogamy scores, that are sub- 
sequently used for distinguishing frustrated regimes of 
quantum spin systems from non-frustrated ones in th ex- 
periment. As we will see below, one of these measures is 
monogamous for pure three-qubit quantum states, while 
the other is not. 

Monogamy Score for Negativity Squared. - Negativ- 
ity is an important quantum correlation measure for 
two-party quantum states [27] . It is defined within 
the entanglement-separability paradigm. The negativ- 
ity, -ZV(i 2 , of an arbitrary bipartite quantum state pi 2 is 
the absolute value of the sum of the negative eigenval- 
ues of the partial transposed state p^J i where the partial 
transposition is taken with respect to Alice (1). The im- 
portance of this measure arises from the fact that if a 
partially transposed bipartite quantum state has nega- 
tive eigenvalues, the state must be entangled [30] . 

Replacing Q in Eq. (JlJ by the squared negativity, we 
obtain the monogamy score for the negativity squared, 
given by 

=^l 2 (23)-^12-^l 2 3- (2) 

Recently, it has been shown that the negativity squared is 
monogamous, i.e. 6^2 > 0, for all three-qubit pure states 
[26] . In this paper, we measure <5/v 2 m both the frustrated 
and non-frustrated regimes. For ease of reference, we 
call the monogamy score for negativity squared as the 
"entanglement monogamy score" . 

Quantum Discord and Monogamy. - Let us now de- 
fine another bipartite measure of quantum correlation, 
and importantly it does not belong to the entanglement- 
separability paradigm. Quantum discord is defined as 
the difference between two classically equivalent formu- 
lations of mutual information, when the systems involved 
are quantum [28], and is given by 

D12 = £>(pi 2 ) = /(P12) - J(pi 2 ), (3) 

where -f(pi 2 ) and J(pi 2 ) are argued to be, respec- 
tively, measures of total and classical correlations of 
Pl2- I{pi2) is defined as S(p{) + S(p 2 ) - 5(pi 2 ), 
where S(o~) = — tr(alog 2 cr) is the von Neumann en- 
tropy of a, pi =tr 2 pi2, and p 2 =tripi 2 . J(pu) = 
S{p 2 ) — <S(p2|i), where the conditional entropy S(p 2 \i) — 
min n i ^2 i PiS(p 2 \n 1 ), with the minimization being per- 
formed over all possible rank-one projection- valued mea- 
surements H] on subsystem 1. Here the output state 
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p 2 |ni = njpi 2 n£/tri2(ll£pi2), and the probability 
Pi =Tr(n^p 12 ). The monogamy score for quantum dis- 
cord (referred later as the "discord monogamy score" ) for 
a tripartite state P123 is given by [29J 



5 D = Di(23) ~ D12 - D 



13- 



(4) 



Unlike the entanglement monogamy score, the discord 
monogamy score, Sp, can be both nonnegative and neg- 
ative, even for pure tripartite states |29) . 



III. FRUSTRATED ISING SPIN SYSTEM 

Frustrated spin sytems have attracted a lot of inter- 
est due to their rich phase diagrams. Moreover, current 
technological developments ensure the possibility of de- 
tecting such phases in the laboratory, which is also one 
of the goals of this paper. Frustration can be observed in 
the system consisting of three quantum spin-1/2 particles 
positioned at the corners of an equilateral traingle, and 
having Ising (nearest-neighbor) interactions. The Hamil- 
tonian of this three-spin transverse Ising model is given 
by 

H = h{al + o 2 x + trl) + J(ala 2 z + 0*0* + a-Jo*), (5) 

where h is the strength of the transverse field, J is the 
coupling strength of the Ising interactions, and |J| 3> 
h. The three quantum spin-1/2 particles are denoted 
as 1,2,3. a l x and a\, for i = 1,2,3, are the Pauli spin 
matrices at site i. When J > 0, i.e, the Ising interactions 
are of anti-ferromagnetic type, the system is frustrated, 
whereas when J < 0, i.e, the Ising interactions are of 
ferromagnetic type, the system is non-frustrated. 

We now describe the adiabatic state preparation 
method used to prepare the ground state of this spin 
system in the laboratory. The quantum adiabatic theo- 
rem states that if a system is initially in the ground state 
and if its Hamiltonian evolves slowly with time, it will 
be found at any later time in the ground state of the 
instantaneous Hamiltonian |31j . The Hamiltonian evolu- 
tion rate is governed by the relation 



Kl;t|«0;i>| 



dt 
f(t) 



< e, 



(6) 



where e is a small number, |0; t) and |1; t) are respectively 
the ground and first excited states of the instantaneous 
Hamiltonian H(t), and g(t) is the energy difference be- 
tween the corresponding energy levels. The system stays 
in the instantaneous ground state of the Hamiltonian 
with a probability (1 — e 2 ) 2 . 

The spin system is initially prepared in the ground 
state of the Hamiltonian h[<j\. + <7 2 + Then the 

system is taken to the frustrated regime by adiabati- 
cally increasing J from to | J m ax\ and similarly to the 
non- frustrated regime by increasing J from to — | J m ax\, 
where | J m ax\ 3> h. The system stays in the ground state 



of the instantaneous Hamiltonian with high probability, 
if J is changed slowly enough so that it satisfies Eq. (|6| . 
The energy level diagram of the spin system is shown 
in the Fig. [T] The energy level of the ground state of 
the system is represented by red curve which is marked 
as Eq. The energy level of the only excited state which 
is relevant in the calculation of the adiabatic evolution 
rate of the Hamiltonian is represented by the other red 
curve which is marked as E\ . The energy levels of all the 
other excited states are shown in blue curves. Though 
there are energy levels in between E and Ei, there are 
no possible transitions from the ground state to these ex- 
cited states, as the transition amplitudes (given by the 
matrix elements similar to the numerator of the left hand 
side in Eq. (|6|)) are zero in these cases. 
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FIG. 1. Energy level diagram. The red curves, which are 
marked as -Eo and Ei , represent respectively the energy levels 
of the ground state and the excited state which is relevant in 
the calculation of the adiabatic evolution rate (Eq. |6|). The 
blue curves represent the energy levels corresponding to the 
other excited states. 



IV. EXPERIMENTAL IMPLEMENTATION 

The spin system chosen for the experimental imple- 
mentation is iodotrifluoroethylene (C2F3I) dissolved in 
acetone-Dg. Here the three 19 F nuclear spins act as 
three qubits. The chemical structure of the molecule, 
the chemical shifts of the three fluorine nuclei, and the 
J-couplings between them are shown in Fig. [2j The ex- 
periments have been carried out at a temperature of 290 
K in 11.7 Tesla magnetic field on a Bruker UltraShield 
AV III 500 MHz NMR spectrometer using a QXI probe. 
The 19 F resonance frequency at this field is 470 MHz. 

In the rotating frame, the NMR Hamiltonian of the 
weakly coupled three-spin system is given by 



"Hnmr 



3 



3 

^ 2 

*<7,=1 



(7) 



where are the chemical shifts of the 19 F nuclear spins, 
Jij are the scalar coupling constants between them. The 
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equilibrium density matrix under high temperature and 
high field approximation, is in a highly mixed state, given 
by 



J eq 



i(7 + C A Peq ), 



(8) 



where £ ~ 10 5 is the purity factor and the deviation 
part of the density matrix [52] 



A Peq oc 7 (al+a 2 z +a 3 z ). 



(9) 



Here 7 is the gyromagnetic ratio of 19 F spin. In liquid 
state room temperature NMR, since the preparation of 
a pure state requires extreme conditions, it is a com- 
mon practice to prepare a pseudo-pure state (PPS) that 
mimics the pure state [33J IS] ■ We have used the spa- 
tial averaging method to prepare the |000) PPS from the 
equilibrium [33]. 
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FIG. 2. Chemical structure of the molecule (left) and the ta- 
ble of Hamiltonian parameters (right). In the table, the diag- 
onal elements are the chemical shifts Vi (in Hz) of the fluorine 
spins and the off-diagonal elements are the scalar coupling 
constants Jij (in Hz) between them. 

In the experimental implementation, the Hamiltonian 
Hit) is discretized into M+l steps as H (jjT), where T 
is the total duration of the adiabatic evolution and m 
goes from to M. The unitary operator corresponding 
to the each step is given by 

U m = cxp(-tH(§T)At) = exp(-i[fc(<7* + a 2 x + u\) 
+J{^T){alal + a 2 z al + ala 3 z )]At), 

(10) 



where At — 



M+l 

(M + l)th step is given by 



The final unitary operator after the 



M 



u=~[[u n 



(11) 



The unitary operator for each step can be approximated 
to second order in At by using the Trotter's formula 

x ex P HJ(§ T){a\al + o 2 z o\ + *la 3 z )At) 
x expH/^ + ul + a 3 x )f) + 0(At 2 ). (12) 

In the experiment, the value of hAt was set to w/21 and 
that of J(T)At to 7r/4 and — it /A in the frustrated and 



non-frustrated cases respectively. The value of jj was 
increased from to 1 in 21 steps, in both the frustrated 
and non- frustrated cases. Considering the energy gap 
between the ground state (Eq) and the relevant excited 
state (Ei), increasing J(t) linearly is not time optimal. 
To achieve a time-efficient adiabatic evolution, we used a 
sine hyperbolic variation in J(t). 

The ground state of the Hamiltonian at time t = 0, i.e. 
of h(al + a 2 x + al), is | ) with |-) = ^(|0) - |1». 

This state was prepared from the |000) PPS by apply- 
ing a 5 rotation with respect to the —y axis on all the 
three spins. This rotation was realized by a numeri- 
cally optimized amplitude and phase modulated radio 
frequency (RF) pulse using GRadient Ascent Pulse En- 
gineering (GRAPE) technique [35]. As evident from Eq. 
(Ill, the k th step unitary operator is a product of k uni- 



tary operators (Uks) acting on the initial state. For effi- 
cient implementation of the adiabatic evolution, we gen- 
erated GRAPE pulses by cascading these unitary oper- 
ators. The length of these pulses ranged between 2 ms 
to 30 ms. All the GRAPE pulses were optimized such 
that they are robust against RF field inhomogeneity and 
the average Hilbert-Schmidt fidelity of all these pulses 
are greater than 0.995. Quantum state tomography of 
the full density matrix was performed after every second 
step in both the frustrated and non-frustrated regimes. 
The reconstruction of the full density matrix involves an 
optimized set of 7 experiments and fitting of the corre- 
sponding real and imaginary parts of the single quantum 
spectra of all the three 19 F nuclear spins [57] . The real 
part of the reconstructed density matrices corresponding 

to the initial state | )( 1 and that of the final 

states corresponding to the last step in both the regimes 
are shown in the Fig. [3] 

To quantitatively evaluate the experimental results, we 
calculate the fidelity (F) of the experimental density ma- 
trices (pexp) with respect to the theoretical density ma- 
trices (p th ), given by 



F 



^(PthPexp) 



2 ) 



(13) 



The fidelity of the initial state | ) ( 1 was found 

to be 0.99 and that of all other final density matrices 
were greater than 0.984. 

The entanglement monogamy score was calculated by 
using the relation in Eq. ([2|. The results are shown 
in Fig. [3^a). It is clear from the figure that the non- 
frustrated regime has higher multipartite quantum cor- 
relations as compared to the frustrated one. Except at 
some extreme points in the non-frustrated regime, the ex- 
perimental results match with the theoretically expected 
curve. To qualitatively analyse the error due to the 
imperfect initial state, we calculated the entanglement 
monogamy score after applying ideal operations on the 
experimental initial state, and the results are also shown 
in Fig. [3]ja). From these results, it is clear that even a 
small loss of fidelity of the initial state can cause large er- 
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FIG. 3. (a) Entanglement monogamy scores. The red-dashed curve corresponds to the theoretically expected results, obtained 
by applying ideal unitary operations on the ideal initial state. The blue circles are obtained by applying the ideal unitary 
operations on the experimental initial state. The blue stars correspond to experimental results, (b) Discord monogamy scores. 
The red-dashed and blue-dash-dotted curves correspond to the theoretically expected results for the pure states and the mixed 
states of the NMR system respectively. These curves were obtained by applying ideal unitary operations on the ideal initial 
state. The blue circles correspond to the results obtained by applying ideal operations on the experimental initial state. The 
experimental results are shown with blue stars. The real part of the experimental density matrices are also shown for the initial 

state | }( 1 (middle one), the final state in the non-frustrated regime (left one), and the final state in the frustrated 

regime (right one). 



rors in the calculation of the monogamy score. The other 
sources of systematic errors in the experimental data in- 
clude decoherence and the non-idealities of the GRAPE 
pulses used to realize the unitary operators correspond- 
ing to the adiabatic evolution. 

We now present our results for the monogamy of quan- 
tum discord. In this case, we make an important change 
in our procedure to calculate the quantum correlation of 
the shared states. Instead of using only the pseudo-pure 
part of the experimental density matrix, as was done for 
calculating the negativity, we used the full mixed state 



density matrix of the NMR system for calculating the 
quantum discord. It is well known that the room temper- 
ature liquid-state NMR systems can have non-zero dis- 
cord, although the purity of these systems is too small 
to exhibit any real entanglement |38j . The discord and 
the discord monogamy score were calculated using Eqs. 
([3]) and Q respectively. The experimental results along 
with the theoretically expected ones are shown in the Fig. 
[3jb). The theoretically expected results are shown for 
both pure and the NMR mixed state density matrices, 
where the latter's pseudo-pure part is the proportional 



() 



to the former. Although the discord monogamy scores of 
the mixed states are very small, their overall behavior is 
very much similar to that of the pure states. The results 
obtained by applying ideal operations on the experimen- 
tal initial state are also shown in Fig. [3^b). Again we 
observe that the experimental data points agree with the 
behavior of the theoretically expected curves. From Fig. 
^b), it is clear that the non- frustrated regime again has 
higher multipartite quantum correlations as compared 
to the frustrated one. For simulating all the theoreti- 
cally expected density matrices, we used the Trotter's 
approximation of Eq. ( |12[ ) and fixed the total num- 
ber of steps for the adiabatic evolution as 21 for each, 
in both the frustrated and non-frustrated regimes. By 
comparing Figs. [3^ a) and j3^b) , we conclude that the 
monogamy scores for negativity squared and quantum 
discord have a similar behavior in both the frustrated 
and non-frustrated regimes. This is despite the fact that 
although the two corresponding bipartite quantum cor- 
relations are defined through widely different approaches 
- one is via an entanglement-separability criterion and 
the other is through an information-theoretic paradigm. 
As noted before, the behavior of each one of the mul- 
tiparty quantum correlation measures is different in the 



frustrated and non-frustrated regimes. 



V. CONCLUSION 

Quantum correlation of separated quantum systems 
are known to be useful in a variety of phenomena. How- 
ever, its detection and quantification remain difficult 
tasks, specially in a multipartite domain. In this work, 
we have been able to use quantum correlations to exper- 
imentally discern between frustrated and non-frustrated 
regimes of a triangular arrangement of quantum spins, 
in a nuclear magnetic resonance system. To attain this 
goal, we have used the behavior of multipartite quan- 
tum correlation measures of the ground states of this sys- 
tem. These multiparty measures are obtained by using 
the concept of monogamy of quantum correlations, which 
puts constraints on the sharability of quantum correla- 
tions. We believe the present study to be important not 
only for understanding quantum correlations in general 
but also in analysing various quantum phase transitions 
in many-body quantum systems, in particular frustrated 
quantum systems. 
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